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We study scattering of three equal mass particles in one dimension. Integrable interactions are 
synonymous with non- diffractive scattering, meaning that the set of incoming momenta for any 
scattering event coincides with the set of outgoing momenta. A system is integrable if the two 
particle scattering matrix obeys the Yang-Baxter equation. Nonintegrable interactions correspond 
to diffractive scattering, where the set of outgoing momenta may take on all values consistent 
with energy and momentum conservation. Such processes play a vital role in the kinetics of one 
dimensional gases, where binary collisions are unable to alter the distribution function. 

When integrability is broken weakly, the result is a small diffractive scattering amplitude. Our 
main result is a simple formula for the diffractive part of the scattering amplitude, when the violation 
of the Yang-Baxter equation is small. Although the derivation is given for 5-function interactions, 
the result depends only on the two-particle scattering matrix, and should therefore also apply to 
finite-range interactions close to integrable. 



I. INTRODUCTION 
A. Background 

In classical mechanics there is an unambiguous notion 
of complete integrability, namely that a system with N 
degrees of freedom must have N (Poisson) commuting 
constants of motion pQ. Finding the appropriate analog 
of this idea in quantum mechanics turns out to be frus- 
tratingly difficult: see Ref. [5] for a recent discussion of 
the pitfalls encountered. 

In Ref. [3], Sutherland persuasively advocates a defi- 
nition of integrability that applies equally well to both 
classical and quantum systems of equal mass particles 
that can move far away from each other (usually in one 
dimension) , out of the range of any interaction between 
particles. In this asymptotic region, the momentum of 
each individual particle is conserved. Scattering pro- 
cesses where the particles begin and end in the asymp- 
totic region are therefore characterized by the initial {ki} 
and final {k^} set of momenta. An integrable system is 
then understood to be one in which these two sets are 
identical. In the classical case, this means that varying 
the 'impact parameters' for the scattering problem does 



not change the set of final momenta, while in the quan- 
tum case, the scattering amplitude vanishes for all other 
assignments of momenta to the outgoing particles. Note 
that in either case, the identity of these two sets does 
not exclude momenta from being exchanged among the 
particles undergoing scattering. In analogy with optics, 
scattering with this character is called non-diffractive. 

Despite being superficially rather different, it is not 
hard to see that these two notions of integrability coin- 
cide in the classical case. Whatever the general form of 
the N constants of the motion for the N particles under- 
going scattering, they must be equivalent to (i.e. func- 
tionally dependent upon) the N momenta in the asymp- 
totic region, which are therefore unchanged after scat- 
tering. In quantum mechanics, the Bethe ansatz is the 
non-diffractive form of the iV-particle wavefunction that 
forms the foundation for the study of integrable systems 
in one dimension. 

For N — 2 particles, energy and momentum conser- 
vation guarantee {ki,k 2 } — {k'l,^}- The distinction 
between integrable and nonintegrable systems therefore 
appears first for N = 3 particles. For an incoming plane 
wave e l ( k i x i+ k 2 x 2+k 3 x 3 ) ^ we wr jt e the scattered wave in 

the asymptotic region as 



^3(xi,x 2 , x 3 ) -t ^ Ap exp [i (knxi + k P2 x 2 + k P3 x 3 )} 
p 
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The first term involves a sum of the TV! = 6 permuta- 
tions of the incoming momenta, with an amplitude Ap 
for each. This part is the Bethe ansatz wavefunction, and 
for an integrable system there is nothing more (strictly 
we must write an expression of this form in each of the 
six asymptotic regions xq\ <C xq 2 <C xq 3i for all six 
permutations Q, but we do not include this extra detail 
for now). The second term, appearing only for nonintc- 
grable systems, is the diffracted wave, a superposition of 
plane waves where the momenta are restricted to a (ID) 
manifold of fixed total momentum P — X)j=i ki and en- 
ergy E = Yh=i |k- Tne amplitude A di g(k[, k' 2 , k' 3 ) is the 
diffraction amplitude. 

For a nonintegrable Hamiltonian the three body prob- 
lem is intractable, even in one dimension. However, a 
system that is close to integrable, that is, whose Hamil- 
tonian deviates only a little from that of an integrable 
model, is expected to display a small amount of diffrac- 
tive scattering. While the primary goal of this work is to 
obtain the diffraction amplitude for three particles in this 
limit, we also seek an understanding of why diffraction is 
sometimes absent. In the theory of integrable systems, a 
distinguished role is played by the Yang-Baxter equation, 
a relation obeyed by the two particle S'-matrix of an in- 
tegrable Hamiltonian. In deriving our result, we will see 
how the violation of the Yang-Baxter equation gives rise 
to diffraction. 

Our result yields valuable insight into how the defin- 
ing characteristic of integrable systems breaks down upon 
changing the Hamiltonian. A more practical motivation 
is provided by recent experiments on one dimensional ul- 
tracold atomic gases [I] , showing essentially no relaxation 
towards an equilibrium state. Recall that in a three di- 
mensional gas, the dominant process of equilibration is 
binary collisions between gas particles, whose effect on 
the evolution of the distribution function is described by 
the Boltzmann equation. In one dimension, such colli- 
sions result in k\ — k[, k 2 = k 2 , or &i = ^2j ^2 = k[, 
which does not alter the distribution of particles in mo- 
mentum space and therefore cannot lead to equilibration. 
Evidently diffractive scattering is required, and three par- 
ticle collisions will be the most important at low density 
[51 15] . We defer to the Conclusion further discussion of 
kinetics due to three particle collisions. 

To what extent are one dimensional atomic gases de- 
scribed by integrable Hamiltonians, or those close to in- 
tegrable? As shown in Ref. [TJ, tight confinement in the 
two transverse directions allows the interaction between 
a pair of atoms to be described by a J-function 

V(xi - x 2 ) = gu8(xi - x 2 ). (2) 

A gas of atomic bosons therefore provides (ignoring any 
confining potential along the length of the gas) a real- 
ization of the Lieb-Liniger model of the ID Bose gas 



[5], soluble by the Bethe ansatz. To destroy integrability 
one may add external potentials, introduce a more com- 
plicated (finite range) interaction between the particles, 
or consider multiple species. While species with differ- 
ent masses certainly lead to diffractive scattering (as the 
above considerations should make clear) , a more relevant 
situation in ultracold physics is to consider different in- 
ternal states of the atoms (e.g. different hyperfine states), 
in which cases all masses remain identical. However, the 
interparticle interaction is dependent upon the species 
involved, and this again leads to a nonintegrable Hamil- 
tonian. Since the variation in the interaction strengths 
is typically on the order of a few percent, this provides 
a natural setting for the question of how weak violations 
of integrability give rise to diffractive scattering. 



B. (5-function potentials and the relation to 
diffraction from a wedge 

Motivated by the above discussion, we take as our prin- 
cipal example the three particle Hamiltonian 

H=p{+pl+pl 

+ 9i25{x\ - x 2 ) + gi3S(xi - x 3 ) + g 23 S(x 2 - x 3 ) (3) 

All masses are equal to 1/2, and we will assume with- 
out loss of generality that all three particles are distin- 
guishable - the scattering amplitude for the case where 
two particles are identical bosons or fermions can be con- 
structed later by symmetrizing or antisymmetrizing the 
solution. The simplicity of Eq. ^ is deceptive: other 
than the case gi 2 = g\ 3 = g 2 3 soluble by the Bethe 
ansatz, only a few other special cases have been solved 
to date [5lU2). The general formulation of the problem 
given in Refs. |12[ 113] is forbiddingly complex, involv- 
ing the solution of a system of functional equations (see 
Eq. (15)). In contrast, we are interested in finding the 
diffraction amplitude when g\ 2 ~ 1713 ~ 1723- Numerical 
treatments of various aspects of the scattering problem 
can be found in Refs. [I3Ml6| . 

For Eq. (3]), the result could be obtained using per- 
turbation theory, starting from the Bethe ansatz form 
of the wavefunction. However, such an approach makes 
no explicit connection to the violation of the Yang- 
Baxter equation, and obscures the fact that the diffrac- 
tion amplitude can be expressed using only two par- 
ticle scattering data. This last feature strongly sug- 
gests that the result is applicable to finite range inter- 
actions close to integrable and not just to ^-function 
interactions. As an example, the interaction potential 
V(xi — x 2 ) — <?i2 sinh~ 2 (a;i — x 2 ) is integrable 3J), so a 
small deviation from this potential will introduce a cor- 
respondingly small amount of diffractive scattering, de- 
scribed by our main result, Eq. (32). 



The solution of the three body problem is greatly fa- 
cilitated by the use of the following picture [5] . By mo- 
mentum conservation the wavefunction has only trivial 



FIG. 1. Geometrical description of three particle scattering in real space. (Left) Particles interact on the three planes defined 
by Xi = Xj. (Right) Projection along the centre of mass motion in the (1, 1, 1) direction. The six sectors correspond to different 
ordering of the three particles on the line, given by a three digit code. 



dependence on the (1, 1, 1) direction in (xi, x 2 , x 3 ) space, 
corresponding to center of mass motion, and we may re- 
strict our attention to the plane perpendicular to this di- 
rection. In this plane, the surfaces x\ = x 2 , x\ = X3, and 
x 2 = S3 appear as three lines at angle 7r/3 to each other, 
cutting the plane into six sectors corresponding to the 6 
possible permutations of the positions of the particles on 
the line (Fig. [T]). The Hamiltonian describing propaga- 
tion in the plane is just the two dimensional Laplacian, 
with the (^-function interactions corresponding to bound- 
ary conditions at the surfaces. We note that the effect 
of different masses may be incorporated into this picture 
by rescaling the spatial coordinates, so that the kinetic 
energy in the centre of mass frame is isotropic, at the 
expense of altering the angles between the three planes. 



C. The Yang— Baxter equation and its geometrical 
meaning 

Let us try to describe the above two dimensional wave 
problem using geometrical optics. We first find the trans- 
mission and reflection coefficients for each of the three 
surfaces, which are defined by the following two particle 
wavefunction 

^ 2 {Xi,Xj) = 



A straightforward calculation using the Hamiltonian 
Eq. (|3| yields 



Uj{ki, k 2 ) 
rij(ki,k 2 ) 



ki - k 2 



k 1 -k 2 + ig i:j 
ki-k 2 + igij ' 



(5) 



Eq. Q describes a collision in which transmission leads 
to the particle with momentum k\ overtaking the parti- 
cle with momentum k 2 . When tracing rays in the two 
dimensional picture, then, this means that k\ — k 2 — 
2fcsina > 0, where < a < it is the angle between the 
ray and the plane, and k = JE — P 2 /3 is the magni- 
tude of the momentum in the centre of mass frame. We 
will therefore write the transmission and reflection as a 
function of a as 



nj(a) -- 



2fcsina 



2k sin a + ig^ 



(6) 



2k sin a 



m 3 



With this convention fixed, we consider different ray 
trajectories with the same arrival and departure angles 
(Fig. |2|. All three rays illustrated have the same path 
length, so the difference in their associated amplitudes 
arises entirely from the reflections and transmissions each 
experiences. For the outgoing wavefront to the right of 
the dotted line, two trajectories contribute, depending on 
whether the wave is first transmitted or reflected, so the 
overall amplitude is the sum 

ti2{a)ri 3 (ir/3 + a)ri 2 {n/'i - a) 

+ r 12 (a)r 23 (ir/3 + a)ti 2 (ir/3 - a), (7) 

with < a < 7r/3. To the left of the dotted line, only 
one trajectory contributes, with amplitude 



r23(7r/3 - a)t 12 (n/3 + a)r 13 (a). 



(8) 



// the amplitudes in Eq. ^ and Eq. ^ are equal, it 
is plausible that the outgoing wave in sector 213 can 
be written as a plane wave with amplitude equal to 
this common value. Returning to the three particle pic- 
ture: if the incoming wave in sector 123 corresponds to 
e i(k 1 x 1 +k 2 x2+k: i x 3 ) ^ ou tgoing wave in sector 213 cor- 
responds to J\. 2 3ie l ( k2Xl+k3X2+klX3 \ where the subscript 
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FIG. 2. Three rays arriving in sector 123 at the same angle, and departing in sector 213, again at the same angle. The paths 
(a) and (b) contribute to the outgoing wave to the right of the dotted line, while to the left only (c) contributes. 



231 on the amplitude indicates how the momenta have 
been permuted. Comparing with Eq. ([I]), we see that 
this corresponds to one term of the Bethe ansatz wave- 
function. 

For the outgoing rays corresponding to the other per- 
mutations of the momenta, one could draw similar sets 
of trajectories. Once again, if equality holds between 
amplitudes for rays contributing to different parts of the 
outgoing wavefront, it seems plausible - and we will show 
explicitly later - that the Bethe ansatz gives the complete 
form of the wavcfunction. That is, there is no diffraction. 

The required equality of Eq. Q and Eq. for ar- 
bitrary a is (one component of) the Yang-Baxter equa- 
tion. Evidently, it is unlikely to be satisfied for an arbi- 
trary set of reflection and transmission coefficients. How- 
ever, for Eq. ^6j, it is satisfied when (and only when) 

9l2 = 9l3 — 923- 

What happens in the general case? In qualitative 
terms, the wavefronts to the right and left of the dot- 
ted line in Fig. [2] will not 'match', having different am- 
plitudes. As we move away from the geometrical optics 
limit, we expect diffraction to smear out this disconti- 
nuity, which sheds light on the connection between the 
violation of the Yang-Baxter equation and the appear- 
ance of diffractive scattering. We next turn to the main 
tool that will be used to make this connection precise. 




FIG. 3. Contour 7 = 7+ U 7- in the a plane for the Sommer- 
feld integral Eq. 1^ [Tg]. 



be seen by writing it as 



7+-7(-?r)+7--7(f) J-f(-Tr) 



D. Sommerfeld integral 

When the Bethe ansatz does not work, we need a more 
general representation of the wavefunction. This is pro- 
vided by the Sommerfeld integral, which provides a rep- 
resentation of the wavefunction 5 , q(t', (f) in sector Q, ex- 
pressed in polar coordinates [17] 



1 

2tt7 



ikr cos a 



A Q (a + (f))da. (9) 



It is straightforward to verify that Eq. ([9| satisfies the 
2D Helmholtz equation [V 2 + fc 2 ] — 0, as long as 
the integrand vanishes at the endpoints. The contour 
7 must be chosen accordingly. Additionally, we require 
that there is only an outgoing diffracted wave. The choice 
shown in Fig. [3] [18 , has the required properties, as may 



The first integral may be evaluated using the residue the- 
orem, while the second and third pass through the saddle 
points of the integrand at 4> = ±7r, and so may be evalu- 
ated at large r to give [T9"j 



^Res-4 Q (aOU «e 



ikr cos (^4>^q — (f^j 
i(fcr+7r/4) 

V Q ((j>), (10) 



where we have defined the diffraction amplitude T>q {jfi) = 
Aq (4> — tt) — Aq (4> + 7r) , and |^q^| gi yc the locations of 

the poles of Aq{o) that are contained within the closed 
contour (after translation by —(f)). Comparison with 
Eq. ([I]) allows us to identify the first term of Eq. ( fl0| with 
the Bethe ansatz contribution, while the second term is 
the diffracted wave. Note that the saddle point at a = 



5 



would give rise to an incoming wave, the reason for its 
exclusion. 

While one is tempted to think of Eq. Q as represent- 
ing a superposition of plane waves with different wavevec- 
tors, with amplitude Aq{o) at angle a + tt [2U], Eq. (flO|) 



makes it clear that the existence of a diffracted wave is 
intimately connected with the absence of periodicity in 
a. Furthermore, as <j> changes, one pole may move out- 
side of the closed contour, while another a distance 2n 
away moves inside. When Aq(ol) is not periodic the re- 
sulting switching of the residues contributing to the first 
term of Eq. ( 10 1 corresponds to crossing the dotted line 
in Fig. [2j where the amplitude within the geometrical op- 
tics approximation changes abruptly. The resulting jump 
in the wave amplitude on crossing this line - not a true 
discontinuity but smeared on the scale of the wavelength 
- is a distinctive feature of the breaking of integrability 
in the far field. 

We will see in Section |IIB| that in this language the 
Yang-Baxter equation appears as a condition for the pe- 
riodicity of Aq{o). 



E. Outline of this paper 

After this lengthy introduction, let us outline the struc- 
ture of the remainder of this paper. In the next section 
we will obtain a system of functional equations obeyed 
by Aq(o) and show how the Yang-Baxter equation cor- 
responds to periodicity of Aq(o), as well as finding the 
explicit form of these functions in the integrable case. 
This is a vital step in the subsequent derivation of our re- 
sult for weak violations of the Yang-Baxter equation. In 
the case of attractive interactions, two and three particle 
bound states can form. In the integrable case, the colli- 
sion of a two particle bound state with another particle 
does not lead to disintegration, even when kinematically 
allowed. Breaking integrability allows this process to oc- 
cur, and we find the amplitude for this process. In our 



Conclusion we discuss the formulation of kinetic theory 
with three particle collisions. 



II. DERIVATION OF THE MAIN RESULT 

A. The system of equations for Aq{o) 

The wavefunctions in the different sectors are subject 
to the boundary conditions of continuity at the line n = 
Xj between two sectors and the condition on the normal 
derivative 



dn 



.'/./I'M', = Xj). 



In polar coordinates this becomes 



1 9* 



fly*. 



(11) 



(12) 



We are going to substitute this into the Sommerfeld in- 
tegral representation Eq. Q , which we rewrite using the 
symmetry of the contour as 

*Q (r, <t>) = i J e~ ikr cos a [Aq (a + <f>) 

-A Q {cf>- a)]da. (13) 



This allows us to use the nullification theorem proved in 
Ref. [TB], which tells us that if an integral of the form 



-ikr cos a 



f(a) da = 0, for /(a) odd (and obeying 



certain technical conditions at -Moo), then f{a) = 0. In 
this way we can convert boundary conditions on ^q(t, 0) 
to conditions on Aq(o). 

We orient each sector in a standard way, measuring 
the angle of our polar coordinates from the bisector of 
the wedge (Fig. [1]) . The above boundary conditions give 
relations between Aq(o) and Aq'(o) in neighboring sec- 
tors, which we write for the sake of definiteness for sectors 
123 and 213 of Fig. [I] 



A 2 iz(a - tt/6) = t 12 {a)Ai 23 {a + tt/6) - r 12 (a)A 2 i 3 (-a - tt/6) 
Ai 23 (ir/G -a) = t 12 (a)A 2 i 3 (-a - tt/6) - r 12 (a)Ai 23 (a + tt/6) 



(14) 



We can express this equation and the five others arising from the other boundaries in a compact form 



/-A 123 (Tr/3-a)\ 
A 213 (a - 7r/3) 
-A 23 i(tt/3 - a) 
A 32 i (a - 7r/3) 
-A 3 i 2 (tt/3 - a) 
V A 132 (a-Tr/3) J 



a 



t/6) 



/ A 123 (a) \ 
-A 213 (-a) 

A 23 i (a) 
-^321 (-a) 
-4312(a) 
V--4i32(-a)/ 



/-A 123 (a - tt/3)\ 
-4213(tt/3 - a) 
-A 23 i (a - 7r/3) 
-4321 (tt/3 - a) 
-A 3 i 2 (a - 7r/3) 
V A 132 (Tr/3-a) J 



S 2 (a- tt/6) 



(A\ 23 (-u 
-A 213 (a 
A 23 i(-a 
-A 32 i(a 
A 3 i 2 (-a 

V--4i32(a 



(15) 



() 



where the matrices Si t2 (a) are 



Si (a) = 



/r 12 (a) t 12 (a) 
*i2(a) r 12 (a) 



V 



S 2 (a) = 













\*23(a) 















^23 (") 
t 23 {a) 









f"i3(a) *i3(a) 

*i3(") ri 3 (a) 









t 23 (a) 

f"23 (a) 













_ , t 12 {a) 
ti 2 (a) r 12 (a] 









h 3 (a[ 
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ri3{a, ,., 
*l3(aO ri 3 (a)J 



t 2 3(a) \ 







(16) 



(17) 



(18) 




FIG. 4. Interpretation of Eq. ( 14 1 in terms of scattering am 



plitudes. Note that angles are measured from the line bisect- 
ing each sector. 



B. Yang— Baxter equation and Aq{o) in the 
integrable case 

We are now in a position to see how our earlier infor- 



mal discussion of the Yang-Baxter equation (Section I C ) 
reappears in this formalism. By keeping track of the ar- 
guments of the Aq{ot) that appear in Eq. (15) and the 
other five relations, we can trace the amplitude corre- 
sponding to each of the rays in Fig. [5] the same rays we 
traced before. We see that in the first two cases the final 
ray amplitude is ^213 (a — 5tt/6), while in the third it is 
-4213 (ct + 77r/6). As we saw in Section IC the combined 
amplitude for the first two diagrams is 

ti 2 {a)r 13 (Tr/3 + a)ri 2 (ir/3 - a) 

+ ri 2 (a)r 23 (Tr/3 + a)i 12 (7r/3 - a), (19) 

while for the third it is 



r 23 {n/3 - a)ti 2 (7r/3 + a)r 13 (a) 



(20) 



Eq. ( 15 ) have a natural interpretation in terms of reflec- 
tion and transmission (Fig. |4| [3T] . They relate Aq (a) 
at an infinite discrete set of a values, which correspond 
to six different rays (Fig. [6]). Each ray corresponds to an 
infinite set of amplitudes with a differing by multiples of 
2tt. 



Eqs. (15) and (16) are written in the reflection diag- 
onal representation, as the entries of our vectors of am- 
plitudes always correspond to the same sectors. In the 
transmission diagonal representation the entries would 
correspond to the rays of Fig. [6] 



Equality of these two amplitudes is then a necessary con- 
dition for ^213(0; — 57r/6) = A 2 i 3 (a + 7w/6) i.e. period- 
icity of Aq{o). 

To determine that these relations are sufficient seems 
daunting at first, as we have to keep track of six sets of 
amplitudes in six different sectors. However, inspection 
of the rays in Fig. [6] shows that scattering follows a well- 
defined order, which provides a great simplification. 

Each boundary between sectors corresponding to a col- 
lision of two particles has three vertices on it, correspond- 
ing to three different scattering events. Scattering pro- 
ceeds in the same order at each boundary. Thus we only 
need to keep track of six amplitudes at a time, which are 
mapped by the S-matrix into 6 others. After scattering 
three times, the rays are in the sector opposite to where 
they started. Scattering three times more will bring them 
back to their original sectors. 

These six scatterings map the amplitudes as follows 



7 




r/3 



7I~/3 — .. <J" 



FIG. 6. All the rays generated by scattering. The color coding 
makes it clear that each of the six directions is present in each 
sector. 



-A213(~Ot) 

-4 23 i(a) 
-^321 (-a) 
-4.312(a) 
\-A 132 (-a)/ 



Si(a-7r/6) 



/-Ai 23 (n/3 - a)\ 
-4 2 i3(a - 7r/3) 
-^23i(tt/3 - a) 
-432i (a - 7r/3) 
--4 3 i 2 (7t/3 - a) 
V A 132 {a-n/3) j 

/ ^ m (a - 4tt/3) \ 
--4 2 i 3 (47r/3-a) 
^ 23 i(a-47r/3) 
-^ 3 2i(47r/3-a) 
A 312 (a-4n/3) 
\-A 132 {4Tr/3-a)J 

Eq. (21) shows that periodicity 



5 2 (ct-7r/2) 



■S 2 (c»-77r/6) 



5i(a-37r/2) 



erioaicity of the Aq{o) is guaranteed if (it's convenient to shift 
S 2 (a — 57r/3)6>i(a — 47r/3) l S>2(a — 7r)5i(a — 
Noting that 51,2 (a) is 27r-periodic, and Si, 2 (a* + tt) = Sf^ («)> fl 



/ A 123 {a -2tt/3) \ 
-^2i 3 (27r/3 - a) 
^23i(a-27r/3) 
-^ 3 2i(27r/3-a 
^3i2(a-27r/3) 
\-^i32(27r/3-a)/ 

/-^i23(57r/3-a)\ 
-4 2 i 3 (a - 57r/3) 
-^23i(57r/3-a) 
-4 32 i(a - 5tt/3) 
--4 3 i 2 (57r/3 - a) 
V ^i 32 (a-57r/3) / 



5i(a-57r/6) 



5 2 (a-ll7r/6) 



/--4.123(tT 

-4 2 i3 (a - 

-^231 (?T 

-4s2i (a - 
--4 3 i2(tt 
V -4i32(a - 

/ -4i 23 (a - 

--4 2 13(27T 

-423i (a - 
-A 321 (2tt 

-4312 (a - 
V-^li32(27r 



-a)\ 
-tt) 

— a) 
-tt) 
-a) 
-tt) / 

2tt) \ 

-a) 

2tt) 

— a) 
2tt) 
~a)J 



a — > a + 7r/6) 
27r/3)5 2 (a - 7r/3)Si(a) = 1 
~ i / vi iv ' > i.w . (i.iivi • i j ' i ■ /' > — •. j ■■}['■ ) i this is CQuivcilGiit to 
Sx(a - 27r/3)5 2 (a - 7r/3)5i(a) = S 2 (a)S 1 (a - n/3)S 2 (a - 2tt/3), 



(21) 



(22) 



(23) 



which is the Yang-Baxter equation. Parts of this equation 
are trivial. For example, one can stay in the same sec- 
tor only by undergoing three reflections, so the diagonal 
elements of the equation are always satisfied. Likewise 
one can only get to the opposite sector with three trans- 
missions. The nontrivial elements are those that connect 
neighboring sectors, and these correspond to the equality 



of Eq. (19 1 and Eq. (20 1 (and the corresponding relations 



for other neighboring sectors) 

Verifying that the matrices given by Eq. ( 16 ) satisfy 
Eq. (23) for the case gi2 — 513 = 323 is now a straight- 



forward (if lengthy) exercise. 

Now that we have a solution of the Yang-Baxter re- 
lations, we still need to find an explicit form for Aq{o). 
From the discussion of Section |ID( we know that the 
Bethe ansatz form arises solely from the pole contribu- 
tions to the Sommerfeld integral. Thus we seek functions 
with the following properties: 

1. Aq{q) are analytic and periodic with period 27r. 

2. Simple poles located at angles corresponding to the 
direction of the incoming wave, and the five other 
angles connecting by scattering (Note that there 
may be other poles off the real a axis coming from 
the scattering matrix, will will play a role when we 
consider attractive interactions). 

3. Residu es of the poles related by the equations of 
Section [HA] (Eq. (p}). 



Why is a relation between the residues enough to guar- 
antee that Eq. ( 15 ) is satisfied for all al We know that 



in the periodic case the Sommerfeld integral is given only 
by its residues. Thus we can invoke the nullification the- 
orem once more to argue that if the residues have been 
chosen correctly, the solution must be correct. 

To take a simple example, consider the impenetrable 
case gij — > 00. Then all reflection amplitudes = — 1. 
Fixing the sector to be 123, we seek a function with a 
pole at a = 0o with unit residue, corresponding to an 



incoming wave at angle 0o + tt. Using Eq. (15), we see 
that there should be poles with residue +1 at a = 0o — 
2n/3 and O — 47r/3, and all angles differing by multiples 
of 2tt, and poles with residue —1 at a = tt/3 — 0o, it — 
00, and 57r/3 — 0o, an d all angles differing by multiples 
of 2tt. 

One way to construct an analytic function with the 
correct properties is via the function p(a) = \ cot(a/2) = 

E°° „ — — having a set of poles with unit residue at 
a = 2im for integer n. Then the function 

^123 ( a ) =p(a-</>o)-p(a-7r/3+0o)+p(a+27r/3-(/>o) 

— p(a ~n + O ) +p(a + in/3 — O ) — p(a + O — 57r/3) 

(24) 

has the correct poles and residues. Simplifying gives 



3 cos 30o 



(25) 



sin 3a — sin 30o ' 
which one may verify has the desired properties. This 
result was originally obtained by Sommerfeld |17j , and it 
should be clear that it works for arbitrary wedge angles 
by the replacement 3 — > 7r/$ (with diffraction occurring 
when 7r/$ is non-integer, so that -Ai mp (a) is not periodic 
in 2tt). 

The extension to the general case soluble by the Bethe 
ansatz should now be clear. The residues are related as to 
one another as implied by Eq. ( 15 ). Denoting the residue 
by H Q (a) 



KesAQ(a) we have, for example 



72-213(00 - tt/3) 
7£23i(7r/3 - O ) 
7?-32i(0o - tt/3) 
^312(tt/3 - 00 ) 

V^132(0O-7T/3)/ 



= Si (0o - 7r/6) 



/ ftl2 3 (0 O ) \ 
7^213 (-00 ) 

7?-23l(0o) 
72-321 (-00 ) 

7^-312 (00 ) 

V72i 32 (-0o)/ 



(26) 

(Note that the minus signs have disappeared!) and simi- 
larly for the S2 equation, after which one constructs the 
27r-periodic analytic function 



"Jl Q ((j)o)p{a - 0o ) + TIq(it/3 - <po)p{a - tt/3 + O ) 
+^q(0 o - 27r/3)p(a + 2vr/3 - O ) + H q {tt - o )p(a - tt + O ) 
< +TZ Q (<Po - 47r/3)p(a + 4tt/3 - O ) + ^q(5tt/3 - o )p(a + 0o - 5vr/3) Q = 123, 231, 312 
72q(-0o)p(Q! + 0o ) + 72q(-tt/3 + (f>o)p{a + tt/3 - O ) 
+^q(-0 o + 2ir/3)p(a - 2tt/3 + O ) + Kq(-k + (f> )p(a + tt - O ) 

. +^q(-0 o + An/3)p{a - 4tt/3 + O ) + K Q {-5ir/3 + <t>o)p{a - 0o + 5tt/3) Q = 213, 321, 132. 



(the B is for Bethe) As we argued above, Eq. (271 must 
satisfy Eq. ( p~5| ) at all a. Additionally, we have verified 
this using Mathematica. 



reproduces the Bethe ansatz solution with no diffracted 
wave. To describe an incoming wave at angle O + tt in 



By construction, the Sommerfeld integral of Eq. (27) 
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the 123 sector we choose 



Diffractive scattering for weak violations of the 
Yang— Baxter equation 



/ ^123(00) \ 






^213(-^o) 







^231 (</>o) 







7^321 (~4>o) 







72-312 (<&)) 







\K 132 (-<t> )J 




W 



(28) 



To deal with the case where two the particles are iden- 
tical bosons or fermions requires an incoming wave of 
appropriate symmetry. 



We now move on to the case where the Yang-Baxter 
equation Eq. ( 23 ) is not satisfied. This problem is much 
more difficult: we must solve the functional relations 
Eq. (151 without the condition of periodicity in a. For 



this reason, only a few special cases have been worked 
out in detail [5HT2"]. 

However, there is a simple result that we can write 
down with little effort in the case that the Yang-Baxter 
equation is nearly satisfied. Recall that the diffracted 
wave has the form 



J(kr+Tr/4) 



(29) 



where Vq (<f) = Aq (4> — it) — Aq (4> + it). This quantity is 
small, because Aq(o) is almost periodic. Using Eq. (21 1, 
we have that 



(-Ai 23 (it - cf))\ 

-4.213(0 - 7r) 
""4231 (it - 4> 

A 32 i(4> - It) 

-^312 (it - 4> 

V ^132(0 -it) J 



= Si(ct>- 2 V3)S 2 - tt/3)Si (4>) 



( A 123 (<f>) \ 

~A 2 x 3 (-()) 
^231 (4>) 

-A 321 (-cf> 
A 312 (4>) 
-An2(-<f>)J 



(30) 



and 



(-A 123 (-it-^)\ 
^213(0 + it) 

~A 23 l (~1t ~ <j>) 
^321 (0 + It) 

-A 3X2 (-Tt - 0) 

V A 13 2(<P + Tt) J 



5 2 (0)5i(^-7r/3)5 2 (^-27r/3) 



/ A m (<l>) \ 

-A 2 l 3 (-& 

A 23 i(4>) 

-A32l(—(f> 
A 3 l 2 ((p) 

\-A 132 (-cl>)J 



(31) 



are almost equal. In the first approximation, then, we can take 



(V 123 (-<j>)\ 
^213(0) 

£ > 23l(-</> 
©321 (<t>) 
2^312(-0 



[Si (0 - 2tt/3)<S 2 (ct>-Tt/Z)S x (0) - S 2 (0)5! (0 - ^/3)5 2 (0 - 2tt/3)] 



( Af 2 M \ 

A% 13 (-<!> 

A^M 

Af 21 (-^ 
Af 12 (<f>) 
\Af 32 (-d>)J 



(32) 



where we have substituted the Bethe ansatz form of 
Aq(a). 



Eq. ( 32 1 is the main result of this paper. It gives a com- 



pact expression for the amplitude of the diffracted wave, 
valid when the quantity in square brackets (which mea- 
sures the violation of the Yang-Baxter equation) is small. 
Substituting the transmission and reflection amplitudes 
Eq. one finds that this quantity is first order in the 
deviation of the interaction constants from the integrable 



point gi2 — <?i3 = g 23 7^ 0, 00. Starting from zero inter- 
action the diffraction amplitude is bilinear in the interac- 
tion constants, consistent with the perturbation theory of 
Ref. [22]. For interaction constants close to infinite, the 
diffraction amplitude is bilinear in 1/p.y. This is to be 
expected since the interaction potential gijS(xi — Xj) may 
be replaced with — (l/g^j) 8"(x\ — x 2 ) for wavefunctions 
that vanish at coincident points. 
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D. Attractive interactions: scattering to and from 
bound states 



With attractive interactions two and three particle 
bound states appear. A two particle bound state ap- 
pears as a surface wave on one of the boundaries between 
sectors. Motion along the boundary corresponds to rela- 
tive motion of the bound pair and unbound particle. A 
surface wave corresponds to a pole in Aq{q) located at 



through the chain of Eq. (21), because S\{if) is evalu- 



complex 



/;(<) 



Bearing in mind Eq. ( 10 ) we have 



tt/6 + if incoming at angle cf> = tt/6 

, j) j — 7r/6 — if incoming at angle cf> = —tt/6 

7tt/6 — if outgoing at angle <j> = 7r/G 

5tt/6 + if outgoing at angle <j) = —tt/6 

where f > 0, and the signs are chosen in order that the 
wave decays as we move away from the boundary. For 
instance, an incoming wave at angle <j) — tt/6 has the 
form 



-ikr cos(7r/6+i(p- 



-ikr cosh <p cos(0— 7r/6) 

v t>~ ^ r s ^ n h *P sin(7r/6- 



(34) 



which describes a wave with wavevector having compo- 
nents k cosh if inwards along the boundary, and ik sinh f 
perpendicular to it. The second factor is nothing but the 
two body bound state wavefunction, which allows the 
identification (f = — arcsinh(<7/2fc). 

In the integrable case inelastic processes in which a 
bound pair forms or disintegrates cannot occur. Break- 
ing integrability leads to a non-zero amplitude for such 
processes. To demonstrate these facts requires that we 
first find the analog of the Bethe solution Eq. ( 27 1 for the 
motion of bound pairs. 

If we start from an incoming surface wave at <j) = tt/6 



in the 123 sector, Eq. (21) maps the amplitude within 



this sector at the following arguments 

tt/6 + if tt/6 — if — > — tt/2 + if 

— » 5tt/6 — if — » — 7tt /6 + if — > 37r/2 — if, (35) 

before we repeat ourselves (due to periodicity of Aq{o)). 
Inspection of the imaginary parts of the waves reveals 
that only the underlined amplitudes correspond to waves 
that decay appropriately at large r. Thus there must be 
no residue at the other values. 

To check that this is so, we start with 



(36) 



corresponding to an incoming bound state of particles 1 
and 2 at the boundary of sectors 123 and 213. We im- 
mediately run into a problem if we try and map forwards 



( n 12 z{Tr/6 + if) \ 






^213(— 7I"/6 - if) 




1 


^23i(tt/6 + if) 







72-321 (-7T/ 6 - i<P) 







72 3 i 2 (t/6 + if) 







\7£i32(— tt/6 - if)) 




w 



ated at the pole corresponding to the two body bound 
state. So instead we map backwards, e.g. 



S 2 (-27r/3-i V ) 



/ n 123 (3Tr/2-if) \ 

K 213 (-3TT/2 + if) 
72-231 (37r/2 - if) 
K 32 i{-3tt/2 + if) 
72-3i2(37r/2 - if) 
\K 132 (-3TT/2 + if)J 
(37) 

In this way one verifies that the three 'forbidden' argu- 
ments of Eq. (35) indeed have zero residue. 



( 72 123 (tt/6 + *</?) \ 
72-2i3(-7r/6 - if) 
72-231 (tt/6 + if) 
72 3 2i(-7r/6 - if) 
72 3 i 2 (7r/6 + if) 
V72i 32 (-7r/6 - if)) 



With the residues in hand, one constructs the Bethe 
function in the usual way (c.f. Eq. (27)). The ab- 
sence of diffraction means that no disintegration of the 
bound state occurs, although non-trivial rearrangements 
are possible. For instance, the outgoing wave at an- 
gle —tt/6 in sector 123, corresponding to the argument 
— 7tt /6-\-if in Eq. (35 ), describes a bound state of particle 
2 and 3, so that the process 



(1,2) 



(2,3) + l, 



is allowed. Note also that reflection of the bound state 
from the single particle does not occur. 

When the Yang-Baxter equation no longer holds, our 
result Eq. (32) gives the leading order diffraction am- 



plitude describing the disintegration of the bound state 
upon collision with the other particle. 



III. CONCLUSION 

The 'quantum Newton's cradle' [3] presumably owes its 
remarkable features to the distribution function of a ID 
gas being unaffected by binary collisions, as described in 
the Introduction. If an experiment of this type were per- 
formed using a gas consisting of more than one species, 
in which the interaction constants were not all strictly 
equal, the three body diffractive scattering described in 
this work will lead to relaxation, which would provide a 
controlled demonstration of the violation of integrability. 
In conclusion, we briefly sketch how our result is incor- 
porated into the kinetic description. 

Ignoring the possibility of coherence between different 
particle species (i.e. off diagonal elements of the density 
matrix), the state of such a gas is described by distribu- 
tion functions fi{k,x,t), where the index i ranges over 
the different particle species, and the Boltzmann equa- 
tion takes the general form 



dt 



[{/,■}]■ 



(38) 



Binary collisions do contribute to the collision integral 
on the right hand side, because two particles of differ- 
ent species may bounce of each other, changing their re- 
spective /i's. However, the two body collision integral 
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will disappear from the equation for the total occupancy 
J2i fit which can only be affected by diffractive scattering 
of three or more particles. 

As an illustration, the three body collision integral de- 
scribing the change in the occupancy of species 1 due to 
collisions with species 2 and 3 takes the form (restoring 
all dimensionful factors) 



dk2dk^ 



E 



tt/6 



dcj>\V Q (<f>,k)\ r < 



2nm J (2tt)- — 
x [/i(fci)/2(fc 2 )/ 3 (fc 3 ) - fi(K)f2(k' 2 )f 3 (k' 3 )} . (39) 



Eq. ( 39 ) is written for a Boltzmann (non degenerate) 



gas for brevity, and only the momentum arguments of 
the distribution functions are shown. To understand the 
relationship of the angle cf> and wavevector k to the mo- 
menta, recall that momentum and energy conservation 
imply 

h + k 2 + k 3 = k\ + fcn + fco = P 

(40) 



This tells us that in three dimensional momentum space, 
the allo wed va lues lie on the intersection of the sphere of 
radius \J1mE and the plane parallel to the (1, 1, 1) direc- 
tion at a distance P/\/3 from the origin (Fig. M). This 



2mE. 



is a circle of radius k — y/2mE — P 2 /3, and the angle 
and sector Q are determined from the angle between the 
points (fci,fc2,&/3) and (k^jk^k^) on this circle. 

The simplest prediction that we can make on this ba- 
sis is that for a two component Bose gas, the relax- 
ation rate of the overall distribution function will be 
oc nin^ + ri2n\ = (n\ + rt2)"-i^2, where ni,2 are the den- 
sities of the two components. 

We are not aware of many studies of the kinetics driven 
by three body collisions. Ref. [25] investigates the effect 
of such collisions on transport phenomena in ID electron 
systems, but stays within the linearized regime. It would 
be interesting to seek self-similar solutions of the three 
particle Boltzmann equation, describing a spatially con- 
stant but non-equilibrium distribution function, similar 
to those found in other circumstances [231. 
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